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Abstract 

Using the recently introduced boundary form factor bootstrap equations, the form 
factors of boundary exponential operators in the sinh-Gordon model are constructed. 
The ultraviolet scaling dimension and the normalization of these operators are checked 
against previously known results. The construction presented in this paper can be 
applied to determine form factors of relevant primary boundary operators in general 
integrable boundary quantum field theories. 

1 Introduction 

The investigation of integrable boundary quantum field theories started with the seminal 
work of Ghoshal and Zamolodchikov [1], who set up the boundary R-matrix bootstrap, 
which makes possible the determination of the refiection matrices and provides complete 
description of the theory on the mass shell. 

For the calculation of correlation functions, matrix elements of local operators between 
asymptotic states have to be computed. In a boundary quantum field theory there are 
two types of operators, the bulk and the boundary operators, where their names indicate 
their localization point. The boundary bootstrap program, namely the boundary form 
factor program for calculating the matrix elements of local boundary operators between 
asymptotic states was initiated in [2]. The validity of form factor solutions was checked in 
the case of the boundary scaling Lee- Yang model calculating the two-point function using 
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a spectral sum and comparing it to the prediction of conformal perturbation theory. In [3] 
the spectrum of independent form factor solutions was compared to the boundary operator 
content of the ultraviolet boundary conformal field theory and a complete agreement was 
found. Further solutions of the boundary form factor axioms were constructed and their 
structure was analyzed for the sinh-Gordon theory at the self-dual point in [4], and for 
the A2 affine Toda field theory in [5]. In the recent paper [6j the validity of the form 
factor solution conjectured for the unique nontrivial boundary primary field in scaling Lee- 
Yang model was tested against truncated conformal space and a spectacular agreement 
was found. 

It is clear from the discussion in [3] that the most interesting open problem of the 
boundary form factor bootstrap is the identification of the operator corresponding to a 
given solution. For example, in sinh-Gordon theory there are infinitely many form factor 
solutions with minimal growth at large rapidities, which can be attributed to the presence 
of exponential boundary fields. The task undertaken in this paper is to map the space 
of such minimal solutions, and make their correspondence with exponential fields more 
precise. It is shown that specific solutions can be selected inside this infinite family, such 
that their scaling dimension agrees with the prediction of conformal field theory, while 
their normalization matches the available results on the vacuum expectation value of ex- 
ponential operators. The construction of these solutions can be generalized to determine 
form factors of relevant boundary primary fields in any model where the form factors of 
relevant operators in the bulk theory are known. 

The outline of the paper is the following. In section 2 the necessary information about 
boundary sinh-Gordon theory is presented. Section 3 describes the construction of the 
form factors solutions which are conjectured to correspond to boundary exponential fields. 
In section 4 their ultraviolet dimension and normalization is evaluated as a series expansion 
in the bulk parameter of sinh-Gordon theory, and is shown to be consistent with known 
results. Section 5 is reserved for the conclusions. 

2 Boundary sinh-Gordon theory 

The sinh-Gordon theory in the bulk is defined by the Lagrangian density 



It can be considered as the analytic continuation of the sine-Gordon model for imaginary 
coupling P = ib. The S-matrix of the model is 



1 m 




B 




where 



sinh \ {0 + inx) 




{x)g{l - x)g 



sinh 6 + i sin irx 



sinh ^ {9 — inx) 



sinh 6 — i sin rrx 
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The minimal bulk two-particle form factor belonging to this S-matrix is [7] 



f{e)=Afexp 



dx 



sm 



X 



XilTT — 



2tx 



sinh ^ sinh(l — -|)| sinh | 
sinh^ X 



(2.1) 



where 



exp 



-4 



dx sinh ^ sinh(l - f ) f sinh ■ 
X sinh^ X 



(2.2) 



It satisfies j(Q^ 5) — > 1 as — > oo, and approaches its asymptotic value exponentially fast. 
Sinh-Gordon theory can be restricted to the negative half-line with the following action 



A 



dt 



dx 



1 TTI 



(2.3) 



+ 



cosh( -($(0,t)-<l>o; 



which maintains integrability [T]. The corresponding reflection factor depends on two 
continuous parameters and can be written as [8] 



B 
1 



B 
1 



E-1 



F-1 



(2.4) 



It can be obtained as the analytic continuation of the flrst breather reflection factor in 
boundary sine-Gordon model which was calculated by Ghoshal in [9]. The relation of the 
bootstrap parameters E and F to the parameters of the Lagrangian is known both from 
a semi-classical calculation [H [10] and also in an exact form in the perturbed boundary 
conformal field theory framework |11) . 



3 Boundary form factors in sinh-Gordon theory 
3.1 The boundary form factor axioms 

The axioms satisfied by the form factors of a local boundary operator were derived in [2] 
and are listed here without much further explanation. Let us assume that the spectrum 
contains a single scalar particle of mass m, which has a two-particle S matrix S{9) (using 
the standard rapidity parametrization) and a one-particle reflection factor R{9) off the 
boundary, satisfying the boundary reflection factor bootstrap conditions of Ghoshal and 
Zamolodchikov pj|. For a local operator 0(t) localized at the boundary (located at a; = 0, 
and parametrized by the time coordinate t) the form factors are defined as 

out{Oi, 6'2, • • • , 6',-„|C(t)|6'i, 62, 



■ ■ ■ 1 "n/in 

id n' n' .on n \ -imt{Y. cosh 9 i-^. cosh B.) 
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for > > • • • > > and 9[ < 62 < . . . < 9'^ < 0, using the asymptotic in/ out 
state formalism introduced in [12]. They can be extended analytically to other values of 
rapidities. With the help of the crossing relations derived in [2] all form factors can be 
expressed in terms of the elementary form factors 

OMt(0|C(0)|6'i, 6^2, ... , On)in = F^iOl, ^2, • • • , On) 

which can be shown to satisfy the following axioms: 
I. Permutation: 



II. Reflection: 

III. Crossing reflection: 

F^{e,, 62,..., e^) = R{nx - e^)F^{2ni - ^1, ^^2, • • • , ^n) 

IV. Kinematical singularity 



-I Res Ff+2(^ + 0', 01, . . . , = - J] s{e - e,)s{e + 0,) j (^i, . . . , 6^, 

V. Boundary kinematical singularity 

n 

Res F^_,,{e + y , 01, . . . , = I (1 - n ^(f - (^1, 



1=1 



where g is the one-particle coupling to the boundary 

There are also axioms for singularities corresponding to bound states (bulk and bound- 
ary), but they are not needed in the sequel. 

The general solution of the above axioms can be written as [2] 

F^e,, 02, ... , 0„) = !^'m"""i!i n ^(^^) n /(^^ - ^^^/(^^ + (3-2) 

i<j 

where the Qn are symmetric polynomials of its variables, and the minimal one-particle 
form factor is given by 

riO)= ...n --''v^-^, •• , j = liE-l) 7' = ^(i^-l) (3.3) 

(smh0 — zsm7)(smh0 — «sm7'j 2 2 



with asymptotics r{6) ~ 1 when 6 ^ oo, wherq^ 

dt 
T 



u{9) = exp 



1 , i 

-—r — 2 cosh - COS 

sinh ^ 2 



sinh 2^ + sinh (l — f) f + sinh | 



sinh t 



and 



4sin7r5/2 



n/2 



(3.4) 



is a convenient normaUzation factor. Using the results of [21 [3] it is easy to derive the 
recursion relations satisfied by the polynomials Qn- 

/C: Q,{-y,y)=0 

Qn+2i-y,y,yi, ...,?/„) = 

(y^ - 4 cos^ 7)(?/^ - 4 cos^ -f')Pniy\yi, . . . , ?/„) Qniyi, ■■■,yn) forn > 
B: gi(o)=o 

Qn+i{0,yi, . . . ,yn) = 

4cos7cos7'5„(?/i, . . .,yn)Qn{yi, ■■■,yn) for n > 



where 



Bniyi,...,yn) = 4g-^^ (l/^-2sin^^ - JJ (^|/i + 2sin 



(3.5) 



(3.6) 



and 



^n(l/bl,---l/n) 



with the notations 



2(y+-l/-) 



,i=l 



i=l 



y+ = ijjz \ ijj ^ z ^ 



(3.7) 



(3.8) 



y- 



U) z + tuz 



with the auxiliary variable z defined as a solution of y = z + z"^ (i.e. writing y = 2 cosh6' 
one has z = e^). 

The two-point functions can be computed from a spectral representation: 



PABimt) = (0|^(t)i3(0)|0) 

oo ^ 



(3.9) 



n=0 



i>...>6»„>0 



d9i . . . dOn /n (6*1, ... , On) exp j — imt^^ cosh 9i 



i=l 



^Notice that the normalization used here differs from that in the earlier papers [2l[3]. 
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where time translation invariance was used and 

/n (^1, • • • , On) = F:^{Oi-, • • • , 9n)^F^{9i, . . . ,9n) 

The fn are symmetric (and also even) in all their variables, therefore the spectral expansion 
(I3.9p can be written in the following form for the Euclidean two-point function 

^ -I -I roo roo i*co 

The operators of interest can be classified according to their scaling dimensions, which 
means that the two-point function must have a power-like short-distance singularity 

PAB{mT) = + . . . (3.10) 
where A^g is an exponent determined by the ultraviolet scaling weights of the local fields. 

3.2 The family of minimal solutions and the cumulant expansion 

In the earlier paper [3] it was shown that there exist an infinite family of solutions for 
which the polynomials Qn have the minimum possible degree 

n{n + l) 
deg Qn = 

These can be thought to correspond to the exponential operators 

^a<S>{t,x=0) 

of which only a countably infinite number is independent since they can all be expanded 
in terms of powers of the field: 



k=l 



For any minimal solution the corresponding form factor has a finite limit when all rapidities 
are taken to infinity simultaneously 



Fniei + A, ^2 + A, . . . , ^„ + A) ^ F„(^i, 02,..., On) + 0(6-^) (3.11) 

which means that every multi-particle contribution in the spectral expansion (13. 9p individ- 
ually behaves as 

with the naive scaling dimension 5 equal to [21 [3]. 
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However, the true scaling dimension in general turns out to be different from due to 
logarithmic corrections in the individual multi-particle contributions, which may sum up 
to give an anomalous dimension. It can be computed using the cumulant expansion of the 
logarithm of the two-point function [I3] (for a very nice discussion see also Consider 
the conformal operator product expansion 

hi 

where h_A and hs are the ultraviolet weights of the fields A and i3, while the hi are the 
weight of the Oi. It is obvious that 

2A_4B = h_A + his — haan 

where /imin is the minimum of the weights hi of the operators Oi appearing in the expan- 
sion. Let us suppose that the limiting function (13.111) satisfies an asymptotic factorization 
property of the form 

Fn{9i, . . . ,9k, 9k+i + \, . . . ,9n + X) = Fk{9i, . . . , 9k)Fn-ki9k+i, . . . , 6'„) + 0(e"^) (3.12) 

both for the form factors of A and B. The leading n = term in the spectral expansions 
is a constant given by the vacuum expectation value of the field, which is assumed to be 
1. In the case of bulk form factors (13.121) with the particular normalization given above 
entails that the operator has a unit vacuum expectation value [I5]. For boundary form 
factors this can be carried over as a reasonable assumption which was used previously to 
normalize the minimal form factor of the boundary Lee- Yang model in [2]. This assumption 
is reasonable from the consistency of the arguments in this paper and is also checked 
explicitly in subsection 4.2. 

Under the above assumptions, the logarithm of the correlation function can be written 

as 

^ 1 -| /»00 /"OO /"OO 

log p^(mr) = V -— — / d9^ d92... rf^„e-"^-S'^°^^^^c„ {9,, ...,9^) (3.13) 

n=l ^ > "'O "'0 Jo 

where the the cumulants of the functions defined recursively by 

/i(^^i)=ci(^^i) , f2{9^.92)=C2{9^,92) + Cr{9^)c^{92) 

M9i, 92, 9s) = cs{9,, 92, 9s) + 01(^1)02(^2, ^3) + 01(^2)02(^1, ^3) + 01(^3)02(^1, ^2) 
+01(^1)01(^2)01(^3) 

Defining 

Cn{9i, ...,9n)= lim o„(^i + A, . . . , ^„ + A) (3.14) 

A— >oo 
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it is easy to see that the functions c„ depend only on the differences of the rapidities. From 
(I3.12P it is easy to obtain the following property of the asymptotic cumulants 

g„(ei,...,4,efc+i + A,...,e„ + A) ~0(e-^) A; = l,...,n-1 (3.15) 

and also note that 

CniOi + \,...,en + \) = Cn{eu • • • , ^n) + 0(6"^) (3.16) 

It can then be shown that 

2A^B = / de2... dOnCn (0, 92,..., On) (3.17) 

n=l ^' ^ ) •> -oo J-oo 

The derivation of this result is a bit more involved than in the bulk case where the transla- 
tional invariance of the form factor in rapidity space can be used. Let us examine the n = 1 
contribution in more detail. The formulae encountered will also be useful in subsection 
4.2. Consider the integral 

oo in 

— Ci(^)e-"""°'^^ 

Due to (I3.16p . Ci{9) approaches its limiting value exponentially fast 

ci(0) = ci + O(e--^) 
The exponential factor has the property 

rriT 

and therefore in the limit of small r 

^ci(^)e-"^-™^^'' {c^{9) - c,) + ^c.^-^roos^e ^g^^g^ 

which can be evaluated as 



27r Jq 271 Jq 2ti 



constant + —KQimr) 
2tx 

Using the asymptotics 

_K'o(mr) ~ — log mr + log 2 — 7e (3.19) 

(where -^e = 0.577215. . . is Euler's constant) the short-distance exponent is given in this 
approximation as 

2A^ = ^ + . . . (3.20) 
Turning to n > 1, the expansion (13.131) can be rewritten as 

^ 111 /'OO /*oo /»oo 

log p^(mr) = ^7^^ / / de^... rf^„e-'"-^«™^'^^»c„ (^i, . . . , 

n=l ^' ^ ' J-oo J-oo J-oo 



using that the functions /„ and consequently also c„ are even in all their arguments. For 
rriT ^ 1, the support of the integrand is concentrated inside the hypercube 

1^,1 <log(2/mr) 

and in the interior (except for a transient shell-like region whose thickness is 0(1) inde- 
pendent of r) 

Q-mTj2, cosh Oi ^ ^ 

Using the properties (I3.15P and (I3.16P c„ can be replaced by c„ up to 0(1) terms. One can 
now repeat the same procedure in the 9i integral as for n = 1, following the line of the 
derivation used in the bulk case [Hj. However, some additional care must be taken since 
the functions c„ are even, and therefore there are 2"~^ independent asymptotic directions 
in which the integrand is non-vanishing: 

c„ (6*1 + A, 6*2 + 0-2 A, ...,9n + o-„A) -> C„ (6*1, Cr26'2, . . . , anOn) 

depending on the choices of the signs cTj = ±1 (the sign of 9i is fixed as only the relative 
signs matter). After a redefinition of integration variables 9i — > atOt this results in 2"~^ 
identical contribution which cancel out n — 1 factors of 1/2. The remaining factor 1/2 can 
be used to map back the 9i integration to the half-line where f l3.19p can be used. 

In the bulk the factor 1/2" are absent, and translational invariance of the form factor of 
a spinless operator implies that c„ = c„. Due to (|3.15l) only the asymptotic region with all 
signs positive contributes, so there is no factor 2"~^ either. Therefore the bulk singularity 
exponent is twice the one in the boundary case (provided that the asymptotic cumulants 
are taken to be identical). However, if the singularity of the bulk two-point function is 
taken to be of the form 

1 

then this factor cancels out in the final formula for A^g and so the bulk cumulant expansion 
is identical to (13.171) . The additional factor of 2 in the exponent is natural since in the 
bulk the conformal weight receives two identical contributions from the left and the right 
movers. 

3.3 The limiting case: bulk form factor solutions in the sinh- 
Gordon model 

Now the task is to resolve the ambiguity in the minimal solution. From [3] it is known to 
originate from the kernel of the recursion relations, which is a homogeneous polynomial of 
order n(n + l)/2 at n-particle level. Let us introduce 

QniVi, ...,yn)= lim A-"("+i)/2Q„(Ayi, . . . , Ay„) 

A—*oo 

which gives the terms in Qn which are exactly of degree n(n + l)/2 (all other terms have 
lower degree). 
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Taking now a solution of the form (13 ■2p , the limiting procedure (I3.1ip gives 
F„(^i,02,...,^n) = lim F„(^i + A,^2 + A,...,^„ + A) 



A— >oo 



n u(^.^^.\ 1 1 / ( - 



i<j 

where Xi = e^* . It was observed in [3] that as a consequence of the boundary form factor 
axioms listed in subsection 3.1, F is always a solution of the bulk form factor axioms 
built upon the same bulk S matrix. For the bulk sinh-Gordon theory, a complete set of 
minimal solutions was found by Koubek and Mussardo in \W\. Let us define the elementary 
symmetric polynomials by 

n n 

l[{x + x,) = J]x"-'a;")(xi,...,x„) 

i=l 1=1 

(t{"^ = if / < or Z > n 

The upper index will be omitted in the sequel, as the number of variables will always be 
clear from the context. Let us also denote 

^ ^ cu" -cu-" sin^ 
\n]- - ^ 



uj — u) sin ^ 

Koubek and Mussardo state that the following family of form factor solutions parametrized 
by the real number k 

. . . , e^) = h /^'^'^'I'Xx^^^ n m - e,) 

(k) 

where the polynomials Pn are given by 

Pi'^ = [k] 

pf) = [A:]detM(")(A;) n>l 

Mif{k) = [i-j + k](72i-j{Xi,X2 ...,Xn) i , j = I, . . . , u - 1 

correspond to the bulk exponential operators 

normalized so that their vacuum expectation value is unity: 

(e'^s*) = 1 
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They support the identification as follows. The above solutions satisfy the asymptotic 
factorization property (I3.12p . derived in [15] for form factors of relevant primary fields. 
In addition, their conformal weights to lowest order in (computed from the two-point 
function) are given by 

Afc = Afc = -— ^ 
87r 

which equal the conformal weights of the exponential operator e'^^* computed in the free 
boson conformal field theory. In section 4 of the present paper more supporting evidence 
is provided for this identification. 

3.4 Resolution of the ambiguity and explicit form factor solutions 

Based on the preceding considerations the ambiguity in the boundary form factor solution 
can be fixed by requiring that 

Q„ = (3.21) 

i.e. 

+ terms of degree strictly less than n{n + l)/2 

For such a form factor the limiting procedure (|3.11l) gives Fn ^ ■ Therefore the limiting 
two-point function cumulants defined in (13.141) will be the same as for the bulk operator 
gfcg* gy comparing the cumulant expansion (13.171) to its bulk counterpart, this means that 
the boundary operator corresponding to this solution has the conformal weight 

eg' 

= 

Sir 

(at least to lowest order in g'^). 

The boundary exponential field e°**^*'°'' has conformal weight 

_^ 
27r 

which leads to the conjecture that the boundary form factors satisfying (I3.2ip correspond 
to the operators 

e 2 

The factor of one-half in the exponent is consistent with the form of the boundary per- 
turbation in the action (12.31) . where a bulk perturbation which is a combination of e^^* is 
accompanied by a boundary perturbation which is a combination of e^^*/^. 

The recurrence relations ( 13. Sp can be solved as follows. The one-particle form factor is 
uniquely given by 

Qf\yi) = Pt^aM = [k]y, 
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The boundary kinematical recursion coefficient (j3.6p can be expressed in terms of tfie 
elementary symmetric polynomials ai in the form 



5„(yi, . . . ,y„) = - 2^ I 2C0S— I cr, 



21 



Supposing that the solution up to — 1 particle level is known, one can write the following 
Ansatz for the n-particle form factor 

Q(f ) = e^e^Bn-iQ^^l^ + a„Pf ) + a^A^f ) (3.22) 

where An is a general linear combination of products of o"/ with total degree strictly less 
than n{n — l)/2, and the first term is understood with the replacement 



and the notation 

ei = 2cos7 , €2 = 2 cos 7' 

was introduced. The Ansatz (13.221) automatically solves the boundary kinematical recur- 
sion B in (I3.5P and therefore the coefficients of An'^ can be determined from the bulk 
kinematical recursion /C. Since the kernel of the recursion is already fixed by (13.221) . the 
solution is unique. It is given by 

Af = 

Af) = [k](el + el - [k]e,e,)a, 

aP = 4sin' ^[k]a^a, + [k]\el + el -[k]e,e,)al(^a, + A sin' ^ 

up to 4-particle level, and it can easily be extended to higher levels using any symbolic 
algebra software. 



4 Evaluating the cumulant expansion 

To provide evidence that the solution indeed corresponds to the exponential operators as 
stated, let us study the short-distance behaviour of the two-point correlators arising from 
it. It was already mentioned that Koubek and Mussardo calculated the ultraviolet scaling 
dimension to order g' by studying the short-distance limit of the two-point function [T6] : 
Babujian and Karowski evaluated it from the cumulant expansion up to two-particle level 
[H]. Here the most singular term in the short-distance operator product of any two of the 
operators in the family of minimal solutions is expanded in powers of g'^ to order g^. In 
addition, the normalization of these operators is also evaluated to order g'. 
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The conformal operator product algebra gives the leading term in the short-distance 
limit of the two-point function as 



with 



where the standard conformal normalization for the exponential operators was assumed. 
Recall that the operator thought to correspond to the form factor solution in the previous 
section is normalized to have unit vacuum expectation value, i.e. it takes the form 



Vk{r) 



1 



kmr) 



:0|e'=t'fM|0) 

so the short-distance expansion can be written as 

where the structure constant C^^'' is given by the ratio of the appropriate vacuum expec- 
tation values. 

The cumulant expansion can be used to evaluate both A^^ and Cj^^'' [lij. The first one 
can be expressed in terms the asymptotic cumulants c^, and as a result the calculation is 
valid both for the boundary and the bulk case, since the asymptotic cumulants Ck are the 
same. This is due to the way the ambiguity in the boundary solution was fixed using the 
bulk solution (see eqn. (I3.2ip ). 

On the other hand, the evaluation of C^^^ involves the full cumulants c^. In the bulk 
theory these are identical to Ck, but for the boundary theory they differ from them. This 
is consistent with the fact that the Ck are independent of the boundary parameters, while 
the vacuum expectation values of boundary fields depend on them in a very nontrivial 
way; therefore the appearance the full cumulant Ck is reassuring. Below it is demonstrated 
that to order the result from the cumulant expansion agrees with the exact vacuum 
expectation values conjectured for the bulk theory in [iTj and for the boundary case in 



4.1 Evaluation of A^; 

Using the cumulant expansion (13.171) . A^; can be evaluated in a systematic expansion in 
powers of g^. The c„ can be obtained as cumulants of the asymptotic spectral functions 

/„^' {e,,...,en)= lim f!:' (^1 + A, . . . , ^„ + A) 

A— >oo 
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which take the explicit form 

\Xi -\- Xj ) 



Xi 6 



From (12 .11] it is obvious that f{9) = 1 + 0{g'^) for any fixed 6' 7^ 0, and the leading term of 

(k) 

the polynomials P„ is given by 

Pf) = fcdetm(")(A;) + 0(/) 
where "mlfik) = {i - j + k)a2i-j{xi,X2 ■ ■ ■ ,Xn) i, j = I, . . . ,n - 1 



while 



As a result 



ckl 2n , ~kl 2n 



In ^ 9 Cn ^ 9 

SO an expansion up to order g"^"- only requires collecting the appropriate terms from cumu- 
lants with up to n particles. The first cumulant can be easily expanded as 



2 " V30 72 2 5671 768 
and using (13.201) gives 

Therefore to recover the intended result it must be shown that all higher order contributions 
cancel, which is demonstrated below up to order g^. 

At order g^ the second asymptotic cumulant is needed. It takes the form 

d2{e,,e2) = \gW{\f{9,-e2f\-i) 

At first this seems to be of order 0(5^^), but some care must be taken. Although it is true 
that f{e) = 1 + 0{g^) for any fixed ^ ^ 0, it also holds that /(O) = 0. Therefore the 
expansion of f{6) in powers of g^ does not converge uniformly in the vicinity of the origin, 
and the series expansion cannot be exchanged with the integration. However, the required 
integral can be evaluated exactly [14]; it has the following small coupling expansioncl 

de {\f{ef\ - 1) = --p2 + o{B^) = 0{g') 

and therefore the two-particle contribution is actually of order 0{g^). 

^In the formula derived in [14] there seems to be an overall sign mistake, which is corrected here. 
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The third cumulant can be written as 



(73 - (Ti(T2 



0-3 - criO-2 



where the are elementary symmetric polynomials of the variables Xj = e^' . The integral 
can be evaluated as 



6 2-K 



27r 



kH + kP kl kl \ 
768 ^ 2567r ~ 3072 J 



which indeed cancels with the term in f l4.2p . 



4.2 Evaluation of C£+' 



C^^*"' is the constant contribution to the logarithm of the two-point function in the short- 
distance limit. It can be expressed as 



g{k,g)g{l,g) 



where Q{k,g) denotes the exact vacuum expectation value computed in |T7] : 



g{k,g) = m''^' {e'^^) 



1 + 



fe2 
2+2b2 



exp 



47? 
dt 



X 



t 



sinh^(2a6t) 



2 sinh^(6t) sinh(t) cosh((l + b^)t) 



and the parameters can be expressed in our notations as 

9 kg 



For small g it can be written as 



and so 



g{k,g) ~ exp 



Clf ~ exp 



2^2 



k'g 



klg\ 



2n 



;-log2 + 7^) 



log 2 + 7^5) 



(4.3) 



On the other hand, it can easily be evaluated from the cumulant expansion in the one- 
particle approximation. For the bulk theory, the first term in (13.181) vanishes exactly and 
sd 

/vr /vr 



-{\og2 -^E)+0{g' 



^ A factor of two was inserted to compensate for the difference between the bulk and boundary rapidity 
integrals. 
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which agrees perfectly with (I4.3p . 

Turning now to the boundary case, the exact vacuum expectation values of boundary 
fields in sine-Gordon theory were conjectured in [18]. With the following action 



A 



bsG 



d X 



^{dtct>f - ^{dAf - 2ficos2f3<f) 
47r 47r 



2fiB / dt COS p{(f){t,0) - (po) (4.4) 



the boundary one-point functions can be written in the following form: 



TTfl 



__oP 

2(1-^2) 



goia,P)gsia,P)gA{a,P) 



(4.5) 



wheri 
ancifl 



exp 



dt 
T 



2sinh^(a/?t) (e 



,(l-/32)t/2 



cosh(t/2) cosh(/?2t/2) - 1 



sinh(/?2t) sinh(t) sinh((l - f3'^)t) 



exp 



/•°° dt sinh^ (a/3t) (2 - cos{2zt) - cos{2zt)) 
T 2 sinh(/52^) sinh(t) sinh((l - (3^)t) 
r /-"^ rft sinh(2a/5t) (cos(2zt) - cos(2^t)) 1 
I Jo Tsinh(/?2t) sinh(t) cosh((l - (3^)t) J 



cosh^vr. = e--^^o^i^ 



cosh nz = e 



2»/3<^o /^l sinvr/j^ 



(4.6) 



provided the operators are normalized as 



^2ia<j>{x) ^—2ia<f>{x) 



gm0(f,O)g-m</.(t',O) 



1 



\x - y\ 
1 



l4a2 



It - 



+ 



+ . . . 



The coupling /i can be related to the mass m of the lightest breather as follows [T9] : 

r \ 1 ^"^'^^ 



7rr(l -/?2) 



2 sin 



T7fP 



2(1-/32) 



2r 



/32 



2(1-/32) 



■^The integral in the exponent of gs contains an additional factor of 1/2 compared to the formula in 
[Is] . The need for such a factor was already noted in [2^, and it was also confirmed in a private discussion 
with ALB. Zamolodchikov. 

^In sine-Gordon theory, z is the complex conjugate of z, but this does not remain true under analytic 
continuation to sinh-Gordon theory. 
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Eqn. fl4.5p can be expanded for small P as 
with 



Due to linearity of gA in k its contribution drops from , so it is not necessary to 
evaluate it explicitly. The sinh-Gordon expectation values can be obtained by taking the 
analytic continuation /3 ig/y/Sii^ under which the operator 

gfcf${T) 

corresponds to putting n = k. Collecting all the terms the final expression is 



klg\ , ^ 1,1 



log 2 + + -TTz coth TTz + -TTz coth TTz) + 0{g ) (4.7) 



47r ' ° ' 2 2 
On the other hand, from the cumulant expansion 



r rif) 

logC^r = ^(log2 - 7..) + / ^ {cm - £i) 



From (13. 3p it is easy to obtain that 



and so 



z sinh ^^(1 — 2 sinh 6^) ^, 

^ ^ = 7-177^ — ■ ■ \( ■ ■ • n + ^ ^ 

(smh — z sm 7) (smh & — i sm 7 j 



de ^ , A;/c/2 /-"^ , r sinh2^(l + sinh2^) 

(ci(0) -ci) = -f^ de> 



Iq 27r 47r (sinh^ 61 + sin^ 7) (sinh^ 6^ + sin^ 7') 

The integral can be evaluated using 

dO TT — 2I7I „ TT TT 

for ~ 77 < 7 < 7: 







sinh2^ + sin2 7 | sin 27] 2 2 
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For definiteness let us suppose that 



TT 



- < 7,7' < 

which means that < -E, F < 1 as a consequence of the parameter identification in (13.31) 
This leads to 



klg^ 1 

^— [log2-7s- -(7r + 7 + 7')cot(7r + 7 + 7') 

-^(7-7') cot(7-7')] 

^[log2 - 7E - -Ti{E + F) cot -7r(E + F) 

-^^^[E-F)coi\'n{E-F)] 



(4.8) 



The expressions (14. 7p and (14.81) match provided 



2*^ - ^> 



(4.9) 



Some care must be taken because there is a branch choice ambiguity in the z parametriza- 
tion which is obvious from (14.61) : in (14. 9p an implicit choice was made. It is only necessary 
to show that there exists one particular choice which is consistent; the correspondence on 
the other branches (i.e. outside the region < i^^, F < 1) can be obtained by analytic 
continuation. 

The consistency condition is that (14. 9p must provide an identification of the reflection 
factor 



R{e) 



1 B 
2+7 



1 



B 



E-l 



(4.10) 



(cf. (12.41) ) of the sinh-Gordon particle to the that of the flrst sine-Gordon breather under 
the analytic continuation [3 — > ig/\/^. The reflection factor of the first breather can be 
written as [9] 



2 2 



1 





'iv 1" 








. ^ 2_ 


e 





(4.11) 



1-/52 

which can be compared with (12.41 ). Under the analytic continuation (3 — > ig/\/^ 

B = -2i 
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and so the factors independent of the boundary condition agree. From the following relation 
between the Lagrangian and bootstrap parameters 



/ sin TT j3'^ 

cos (/J^r^) cosh = /^b W cos (/30o) 



/i 



sin i^P'^rj) sinh = /is W sin (/?0o) 

derived by Al.B. Zamolodchikov pLlJ (see also 120]), comparing to (14.61) the following relation 
is obtained: 

z=—{r]-i§) , z = —{r] + i^) (4.12) 

TT TT 

where an appropriate choice of branches is understood. To the leading orderly in (14. 9p 
and (I4.12P are equivalent to 

TT ' TT 

which makes the reflection factors (14.101) and (14. lip identical. 



5 Discussion and outlook 

With the method developed in this paper, it proved possible to organize the infinitely 
many minimal boundary form factor solutions of sinh-Gordon theory encountered in [3] 
into well-defined families by restricting the highest-growing term to be equal to the bulk 
form factor of a given exponential operator. It was shown that such a solution has the 
correct scaling dimension to order g^, and the claim that its vacuum expectation value is 
normalized to 1 matches with the exact vacuum expectation values conjectured in [18] to 
order g"^. 

The restriction leading to unique determination of the solution was made possible by 
the observation made in [3] that the dominant asymptotic contribution of the boundary 
form factor always solves the bulk form factor equations. It is very important to note that 
the construction always guarantees that the cumulant expansion for the scaling dimension 
of the boundary operator coincides with that of the bulk one which was used as an input 
for the solution of the boundary form factor bootstrap axioms. 

This observation leads to a proposal for constructing form factors of relevant boundary 
operators. In the bulk, form factors of spinless relevant fields can be selected from the 
space of solutions to the form factor axioms by the asymptotic factorization condition 
(I3.12p which was derived in [15]. Using these solutions as an input one can then construct 
solutions for relevant boundary fields whose scaling dimension coincides with the right (or 
equally, the left) conformal dimension of the corresponding bulk field. In fact, the form 

^ Since eqn. (|4.9|) itself is only valid to leading order in it doesn't make sense to take the comparison 
beyond this approximation. 
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factors of the only primary boundary field in the Lee- Yang model computed in [21 [6] could 
also have been derived this way. 

In principle, the agreement of the most singular term calculated from the form factor 
expansion with the quantum field expectation for some operator O leaves open the possibil- 
ity of mixing with other fields of smaller weight (i.e. which are more relevant) in the sense 
that it allows that the form factor corresponds to O with an addition of an unspecified 
linear combination of operators with smaller weight, i.e. to an operator of the form 

where the Oi are operators with scaling dimension less than that of O. This problem 
was already noted in [^. However, for the most relevant boundary field this issue is 
definitely not present which is reassuring for the determination of its form factors using 
the construction proposed above. 

In the sinh-Gordon case the spectrum of conformal dimensions is unbounded from 
below, but the fact that the short-distance limit of the two-point function constructed 
from form factors is consistent with the leading part expected from the conformal operator 
product algebra and the exact vacuum expectation values eliminates some possibilities of 
such mixing. However, it still allows e.g. for 

l>k 

(where it was assumed that k > 0). 

Another restriction comes from the continuation to sine-Gordon theory. The parameter 
k of the solution which gives the corresponding boundary exponential operator as 

gfcf1>(f,0) 

can take any real value, but the operators with k E Z form a closed operator algebra, 
which is the minimal one consistent with the presence of the boundary potential in (12.31) . 
Continuing analytically to the sine-Gordon theory with action 04.41) . the n-particle form 
factors of these operators become the form factors of 

with a many-particle state that consists exclusively of n first breathers. The operators with 
k = ±1 are the most relevant ones in this family, and therefore they cannot mix with other 
operators at all; this property is expected to hold after continuing back to sinh-Gordon 
theory as well. In fact, these would be the operators for which there is the widest spectrum 
of subleading operators to choose from in (15. ip . therefore the absence of mixing in their 
case is particularly striking. 

The above arguments cannot be considered complete, but taken together they make it 
likely that there is eventually no operator mixing for the exponential fields in the boundary 
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sinh-Gordon model at all. The basis provided by the solutions constructed in this paper 
is the only one for which the following property holds: for any member of the basis, the 
limiting form factor (in the sense of (13.111) ) is not only a solution of the bulk form factor 
axioms, but also satisfies the asymptotic factorization condition (13.121) . Note that this 
factorization property is needed for the derivation of the cumulant expansion (13.171) and 
it also ensures that the integrals in the expansion are convergent (cf. eqns. (I3.15P and 
(I3.16P ). Certainly some deeper understanding of the boundary form factor bootstrap is 
needed in order to relate this condition to properties of relevant boundary primary fields, 
and to make the correspondence between form factor solutions and operators more precise. 
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